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ABSTRACT 


A computer program is utilized to examine the effect 
on sound propagation cf the oceanic polar front in the 
Denmark Strait and aiong the Southeastern coast of Green- 
land. Ray traces are computed with the source in both 
the cold and warm current areas for surface thermal 
gradients of 0.036 C/KM to 1.68 C/KM. These gradients 
are created by oe horizontal distances between actual 
oceanographic stations. The results indicate greater 
trapping angles with the source in cold water. A relation- 
ship is reported between the surface thermal gradient and 
the change in trapping angle difference for pairs of runs. 

The location of the oceanic polar front and the effects 
of environmental influences on its strength and seasonal 


movement are described, 
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I. INTRODUCTION 


A. Statement of the Problem 

The term "oceanic polar front" has been used by 
Norwegian oceanographers for over a half century; however, 
the term "front" to describe oceanic boundaries came into 
use in the United States only within the last two decades, 
Fuglister (1954) referred to a frontal zone as separating 
relatively permanent adjacent currents and Cromwell and 
Reid (1956) attempted to formalize a definition for small 
scale thermal fronts. As late as 1961, LaFond proposed a 
more general definition of a front as "The leading edge 
of a border separating unlike water masses...." The 
definition reported by Dietrich (1964) and used during the 
Polar Front Survey conducted during the International 
Geophysical Year 1957-1958 will be utilized here. This 
specifies an oceanic polar front as the boundary separating 
cold, sub-polar water of low salinity from a water mass 
composed of warm, highly saline water originating in sub- 
tropical regions. 

Other than the west wall of the Gulf Stream off the 
east coast of the United States, very few of the oceanic 
frontal zones have been thoroughly examined. Even the 
"Gulf Stream wall" has not been subjected to the precise 
synoptic sampling required and there is still much we do 


not know about its structure and dynamics. 


As may be inferred from the definition above, a 
frontal zone is an area of very complex temperature and 
Salinity structure. The large horizontal gradients found 
in these regions greatly affect the behavior of sound 
waves in both the vertical and horizontal planes, and 
oceanic "abnormalities" such as temperature inversions 
and sub-thermocline ducts frequently are found at the 
fronts. 

Since the oceanic fronts are convergence zones between 
two currents with different physical, chemical and biolo- 
gical properties, changes in water color, wave action and 
other physical parameters often may be observed at the 
boundaries from either aircraft or surface vessels. 

Abundant marine life often is found in the frontal 
regions and the zonal boundaries usually will mark the limit 
of distribution of certain species. The production of 
phytoplankton is usually high at the fronts. Furthermore, 
there is also an accumulation of zooplankton, partly brought 
about by the transport of organisms by the converging 
currents. As a result, frontal areas provide a large 
percentage of the feeding areas of commercially important 
fish species in the northern seas (Shpaykher, 1964). 

Due to the biologic activity and the complex water 
mass structure, the effectiveness of hull-mounted sonars 
may be degraded in these regions and false targets will 


abound, Thus the importance of the fronts is twofold: 


10 


fret secie erfect of fronts relative. to commercial 
fisneries; and second, the effect of increased diffi- 
culty of detection in the military Anti-Submarine 
Warfare problem. 

The purpose of this study is to examine the hydro- 
graphy of the water masses comprising the oceanic polar 
front in the Denmark Strait and to determine frontal 
effects on sound propagation across the frontal axis. 

The hydrography of the area as reported in the liter- 
ature has been compared with seven oceanographic sections 
taken across the polar front in 1955 and 1958. 

In order to examine the frontal effects on sound 
propagation, an assumption was made for each section that 
the two stations lying across the front represented the 
boundary of a different water mass. Sonar ray trace 
programs were then run for each section using assumed 
horizontal separations of ten, six and four kilometers 
between the stations. The bottom was assumed level at 
400 meters in four of the sections, level at 230 meters 
in two of the sections and level at 300 meters for the 
final section, The ray traces were computed using a 
source depth of 25 feet simulating hull-mounted sonar. 
Two sets of traces were computed for each section. In 
the first instance, the source was placed on the warmer 
side and ray traces were computed looking into cold water. 
For the second set, the source was then moved to cold 


water and runs made toward the warmer profile. 
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The ray traces were completed using a two-dimensional 
ray trace program on ths Control Data Corporation Model 
6500 computer at Fleet Numerical Weather Central, (FNWC), 
Monterey (Ayres, et al. 1966), and are described in detail 
in Section iV. For eacn individual station used in the 
ray trace program, temperature and salinity profiles are 
shown, Selected ray traces and their sound velocity 
profiles are also illustrated (Figures 12 to 15). 

Bo The Data 

Since the problem deals with an ocean area con- 
taining severe horizontal temperature snd salinity gra- 
dients, the technique of long-term data averaging did not 
appear to offer a solution. As an alternative approach, 
it was decided to employ quasiesynoptic sections taken 
across the frontal zone in an attempt to describe the 
structure of the front on a seasonal basis, The sections 
utilized were obtained from FNWC, Monterey, and comprise 
those contained in the FNWC Ocean History Information 


Retrieval Systems (FOHIRS) from the ocean area delineated 


in Figure 1, for the years 1955 to the present (Griswold, 
et al. 1968). Only seven ugable sections from the area of 


| 


terest were contained in the data bank. Even though the 


ft 
> 


sections, as snown in Fisure 2, provide good coverage of 
the Strait, it wovid be desirable to have seasonally 


repetitive data for the same location, 





Figure 1. Area of data search 





Figure 2. Location of the sections 


3.3 


The seven sections obtained from FNWC, Monterey, 


provide seasonal coverage as follows; 


Winter January-March a 
Spring April-June Ly 
Summer July-September 2 
Fall October-December 0 


Five of the sections were taken in 1955 and the remaining 
two in 1958. For each section used, temperature and 
salinity profiles are shown in Figure 3. 

The current structure at the front in 1955 was charac- 
terized by above-average development of the East Greenland 
Current and a more weakly-developed Irminger Current 
(Hermann, 1959). Conversely, in 1958 the Irminger Current 
showed its greatest intensification in many years while the 
East Greenland Current was weak (Hermann, 1960). 

These observations are confirmed by Smed (1965) who 
Shows that the yearly anomalies of the surface temperature 
in the area reached their highest value in fifteen years 
during 1958 while 1955 was the second coldest year of the 
decade, 1950-1959. 

Due to the varying position of the sections, seasonal 
influence, and frontal meander, no direct correlation with 
annual variations is possible from the sections. For the 
Same reasons, seasonal correlation was not feasible although 
the two September sections, C and E, do show the lowest 


surface salinities. 
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A problem encountered in attempting to deal with a 
narrow frontal zone concerned the horizontal distance be- 
tween the individual stations which comprise each section. 
These ranges average ten to fifteen miles and provide, at 
best, a rough description of a zone in which surface tem- 
peratures may differ by as much as 6 C across a distance 
of only three nautical miles (Krugler, 1952). 

An examination of this area using closely spaced 
Experimental Bathythermograph (XBT) drops or a thermistor 
chain would be desirable but unfortunately has not been 
accomplished. Mazeika (1968), in his work across the Gulf 
Stream, has demonstrated that excellent approximations can 


be made using non-continuous data. 
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II. THE CHARACTERISTICS OF THE WATER WASSES 

Ae TOC CUT Ieee 

The Bast Greenland Current flows in a southerly 

direction on or near the continental shelf of the eastern 
coast of Greenland (Figure 4), It may be traced from the 
straits between Northern pcenegueaS West Spitzbergen 
as far south as Cape Farewell, a distance of approximately 
3000 km. It is composed, at origin, largely of water 
from the Arctic Basin and may have a characteristic temn- 
perature as low as -1.72 C with salinities from 31.00 0/oo 
to 33.5 ofoo depending largely on the influence of mixing 
and ice thaw or freeze. As the current moves to the south, 
it forms a narrow, shallow band of cold water lying against 
the western edge of the warmer and more saline Greenland 
Sea. According to Chaplygin (1959), the upper boundary of 
the East Greenland Current is always within 55 meters of 
the surface and at no time is it thitkerm than 21 0emevere. 

To the north of fhe Greenland~Iceland Ridge, the East 
Greenland Current is underlain by what several authors refer 
to as the Return Atlantic Current (Figure 5). This cumpent 
forms the western side of the counterclockwise gyre in the 
Norwegian/Greenland Sea area and is usually found at depths 
between 200 and 300 meters with a pense ratte near, 2 C and 
a salinity of 34.75 = 34.90 o/foo (Lee, 1963). 

According to Chaplygin, it is a slow current flowing at 


only 1-2 em ae a, while Sverdrup (1942) gives a figure of 
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Figure 5. Cumrent profile nerth of the 
Greenland-Icelandic Ridge 
(after Gladfelter, 1964) 


be 


25=355 om ae for the core of the more rapid East Green- 
land Current. 

A major portion of the slow Return Atlantic Current, 
or undercurrent, is turned toward the East in an area to 
the north of the Greenland-Iceland Ridge and flows along 
the north coast of Iceland. South of the Greenland-Iceland 
Ridge, the East Greenland Current is bordered and under- 
lain by the Irminger Current water as it proceeds toward 
Cape Farewell. 

To further complicate the thermohaline structure 
near the bottom, there is evidence of intermittent over- 
flow of cold deep water across the Greenland-Iceland Ridge. 
This flow, formed by winter cooling in the Norwegian Sea 
(Harvey, 1961), is defined by Cooper (1955a) as "a series of 
self-contained, cold, heavy balls or boluses of water", 
Harvey (1961) reports a decrease in temperature at the bottom 
of the Denmark Strait in August, 1960 from 2.5 C to 0.1 C 
with a sampling interval of 20 hours. The measured bottom 
current increased from 0.3 kts to a range of 0.54 kts to 1.4 
kts, flowing in a southwesterly direction roughly parallel 
to the trough through the Greenland-Iceland Ridge. 

The Irminger Current is a branch of the North Atlantic 
Drift Current, a part of the Gulf Stream System, and brings 
Warm saline water into the Irminger Sea. As the Irminger 
Current spproaches the southwest corner of Iceland, it 
splits into two branches. The major portion of the water 


moves to the southwest and flows alongside and in the same 


Se, 


direction as the East Greenland Current while a much 
smaller amount courses around the Northern coast of Iceland 
toward the Faroe-[celand Ridge. The Irminger Current water 
may vary seasonally, but many authors assign a temperature 
of 8-10 C and a salinity of 35.20 o/foo (See for example 
Cooper, 1955b) as average values. - This water also has been 
referred to as North Eastern Atlantic Water (Holzkamm, 

et al. 1964). 

Tne Irminger Current is driven by the wind system of 
the Icelandic atmospheric low as will be shown later, This 
low has large year-to-year-variations in its mean intensity. 
Thus, relatively large changes may be expected in the volume 
of flow of the Irminger Current as well as in thé posi#woen 


of its boundaries. 


Be The Vertical Structure 
As may be seen in the section profiles (Figure 3), 

tne vertical gradients are usually very gentle in the warm 
region. Upon entering the colder zone however, the vertical 
sradients of both temperature and salinity increase sharply 
due to the shallow extent of the East Greenland Current 
and the mass of warmer, more saline water wnich underlies 
19% 

Or note is the’ frequency, of sSowmmpncnee of sub-thermo- 
cline ducts found at stations on both sides of the front. 


A preferred depth for these ducts appears to be 75 to 100 


34 


meters; their presence suggests that Variable Depth Sonar 


might be used with excellent efficiency in this region. 
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III. INFLUENCES ON THE STRENGTH AND LOCATION OF THE FRONT 
The strength of the horizontal gradients in the frontal 
area is a function of several interacting environmental 
processes, The major parameters considered in this study 
are: 
a) meteroiogical influences (wind currents at the 
surface causing oceanic current convergence), 
b) ice (affecting the convergence of both wind and 
thermohaline circulation) and 
c) pevren popoeRaene 
These parameters, with the exception of c), vary daily as 


well as seasonally or annually. 


A. The Meterological Influence 
As indicated in the Local Area Forecasters Hand- 
book for Keflavik, Iceland (1963), the area south of 
Greenland and southwest of Iceland is a junction for three 
diverse air masses: 
1) the warm maritime air masses moving northward 
from the central North Atlantic, 
2) arctic air from the arctic Ocean or the Green- 
land ice cap and 
3) continental air masses moving eastward from 
North America. 
These air masses frequently interact to form an atmospheric 
polar front lying on or néar the axis of the oceanic solae 


front. Figure 6 shows the mean seasonal positions of the 
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two fronts. The atmospheric polar front exhibits large 
fluctuations with time and the Location shown is subjective. 
It is difficult to state a causal relationship between 
the fronts, but one may conclude that the surface currents 
which make up the oceanic polar front are greatly affected 
by the wind system. By extension then, we may postulate 
that the position of the oceanic polar front is directly 
related to that of the atmospheric polar front and that it 


reacts seasonally to atmospheric changes. 


B. Ice and Its Influence 

The seasonal distribution of ice along the coasts 
of Greenland and through the Denmark Strait is shown in 
Figure 7 (Climatological and Oceanographic Atlas for 
Mariners, 1959). The ice distribution may vary greatly 
from year to year and even day to day. Meyer (1965) has 
shown that during light ice years, the atmospheric high 
pressure system over Greenland is highly developed and the 
Icelandic low shifted to the south of its long term mean 
position. Conversely, heavy ice years show the [Icelandic 
low shifted to the north. According to Meyer, the ice field 
moves onshore and becomes closely packed as a result of 
northerly, northeasterly, and easterly winds and will be 
pushed well off the coast as a result of westerlies. Figure 
8 from Meyer shows an example of the rapid changes possible. 
On 21 April 1960, the ice edge was 56 miles from the coast. 


Eleven days later, after a period of northeasterly and 
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gure 7 a. Average ice conditions in May 
fter Climatological and Cceanographic 


Atlas for Mariners, 1959) 


—~'t 
| 





Figure 7? bd, Average ice conditions in September 
~wewe-s-> Mean ice concentration 0.8-1.0 

—-—--- Minimum extent, concentrationo,5 or greater 
scuscexe Maximum extent, concentration®.5 or greater 
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2-9 May 1960 
19-21 April 1960 
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36 
Figure 8. Rapidly changing ice boundaries 


off the S.E. coast of Greenland 
(after Meyer, 1965) 
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easterly winds, the ice edge was only 4-6 miles from the 
coast. Meyer further states that the ice edge is rarely 
driven beyond the border of the continental shelf. It 

will be shown in the next section that the oceanic polar 
front normally lies along the edge of the shelf. Since 
only the most massive ice floes and large bergs could exist 
for any length of time in the warm Irminger Current, it 
would apueor the ice distribution strengthens the hypo- 


thesis of front location along the shelf edge. 


C. The Bottom and Its Influence 
The bottom topography of the Denmark Strait and 

Southeastern Greenland is varied. The Greenland-Iceland 
Ridge connects the shelves of Greenland and Iceland at 
depths of 200-300 m, but is cut along its NE-SW axis by a 
600 M channel. The shelf off SE Greenland is twenty to 
fifty miles wide and is heavily glaciated with a discern- 
ible line of terminal moraines found near the break in slope. 
The slope is precipitous, falling from depths of 400 to 
LOOO meters in four or five miles, and then more gradually 
deepening to 3000 meters in the Icelandic and Labrador 
Basins (See Figure 9). 

Along the Greenland-Iceland Rise, the oceanic front 
is more diffuse than on either side of the Rise. It is 
believed by the author that this is due partly to the more 


easterly position of the Irminger Current in this region 
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Bottom morphology of the North Atlantic 


Figure 9. 


1964 ) 


(after Dietrich, 


partly to the topographic deflection of the majority of the 
Return Atlantic Current, 

South of the Greenland-Iceland Ridge, the center of 
the frontal zone tends to lie along the 1000 m contour, 
This is demonstrated by Dietrich and Gieskes (1966) and 
Dietrich (1958). As presented in cross-section by Hermann 
(1959, 1960) and reproduced in Figure 10, the front is seen 
to lie directly over the 1000 m contour. The close prox- 
imity of the 400 m contour to the edge of the slope and 
the surface thermal gradient shown in Figure 11 indicate 
that the shoreward edge of the front closely follows the 
400 meter level along the south-eastern coast of Greenland. 
The depths assumed for the ray trace programs were based on 


this postulation. 
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Figure ll. Temperature along the 
S.=. Greenland Slope, 1966 
(after Dietrich and Geiskes, 1968) 
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IV. FRONTAL EFFECTS ON SOUND PROPAGATION 
A. The Computer Program 

As indieated in» part I.B., sonar ray traces were 
utilized to investigate frontal effects on sound propa- 
gation. The ray trace method was employed since other 
approaches, such as the normal mode technique, would cause 
insurmountable difficulties in this highly complex region. 
No precedents for this study have been found in the litera- 
eure. 

The ray traces were computed using the high precision 
computer program developed at FNWC (Ayres, et al., 1966). 
This program computes the ray position every 1/128th of a 
second of travel or less. AS a ray approaches the surface 
or bottom, the time step is reduced so that the actual 
contact point is computed, Additionally, the program 
corrects for earth curvature (Wolff, et al.). Inputs 
mequirea for each run are: 

a) temperature and salinity at selected depths, 

surface to botton, 

b) range of each profile from sound origin, 

ec) bettonm configuration, 

ad) source depth and 

e) angles of propagation. 

Input depths were selected to provide the computer 


with information on changes in temperature and salinity 


nil 


gradients and tend to be clustered about points of maximum 
change. 

Output data permitted location of the ray in range and 
depth from the source at 60, 150, 300 and 500 feet from the 
surface and whenever the ray touched the surface or bottom 
or executed a maximum depth or eR depth turn. These 


data are tabulated for each set of profiles in Tables 1-7. 


Bo Assumptions 
The principal assumption made in this study of 

the effect of the oceanic polar front on sound propagation 
is that the stations lying on opposite sides of the front 
properly represent the water masses found there and that 
by compressing the profile spacing we are duplicatimg 
thermohaline gradient conditions which may be found in the 
environment but for which no data presently exist. The 
justification for this compression is the high surface 
thermal gradient reported by Krugler (1952), Dietrich, 
et al., (1961), and others. These gradients were not found 
in data available to the author, but do exist. There is no 
technique for precisely locating the front on any secwieg 
and the thermohaline structure present at any station 
justifiably may be assumed to continue toward the boundary 
for some distance, thus increasing the gradients and reduc- 
ingethe compression wmequired, 

Compression of the stations ranged from 1.4/1 for the 


10 KM run for Section C to a maximum of 8/1 for the 4 KM 
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run of Section A. Compressed surface temperature gradients 
ranged from .036 C/KM in the 4 kM run of Section F. The 
maximum gradient of 1.68 C/KM exceeds the surface gradient 
reported by Krugler (1952) of 6 C in 3 NM (which converts 
to approximately 1 C/KM), but is much smaller than the 
gradient reported by Dietrich et al., (1961) of "7 C over 
a distance of a few hundred metres". 

As previously stated, the bottom was assumed level for 
each of the sections, but a depth representative of its 
location was assigned to each. This was done to permit 
comparison of frontal effects from section to section to be 


made without the variable influence of bottom topograpny. 


C. Results 
The compilation of results in Tables 1-7 presents 

the section designator, date, season, and station numbers 
used in the ray trace programs with their locations; 
charted depths and assumed depths. The actual horizontal 
distance between profiles is given, as well as the difference 
in measured surface temperature, AT (SURFACE). Assumed 
separations of 10, 6 and 4 KM were divided into AAT (SURFACE) 
to obtain the DT/DX gradient information. 

The angle trapped is defined as the angular distance 
in degrees between the horizontal or zero degree ray and 
the angle which first touches bottom within the range of 


consideration. 
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For each depression angle (@), 0,1, 2, 3, 4, 6, 8, 
and 10 degrees, the number of surface bounces (x), number 
of maximum depth turns (y) and number of minimum depth 
turns (z) are recorded in the columns below the corres- 
ponding letter, 

The final tabulation gives the maximum depth attained 
by the ray (in feet) and the horizontal range from the 
source (in kiloyards) at which this depth is reached. 

Upon examination of these data in Tables 1-7, several 
points of interest became apparent: 

1) In all runs where the assumed surface thermal 
gradient was greater than 0.20 C/KM (36 runs) the cold-to- 
Warm propagation tended to produce more surface bounces 
than the warm-to-cold runs. The ray density is greater in 
the near-surface layers in the cold-to-warm runs. MThis is 
further demonstrated by the maximum depths attained by the 
rays, in that the rays from cold-to-warm were more shallow 
than the corresponding warm-to-cold ray. This comparison 
is illustrated in Figure 1lé, 

2) Conversely, in those cases where the assumed 
surface thermal gradient was less than 0.2 C/KM, the warm- 
to-cold propagation produced more surface bounces and the 
cold-to-warm rays reached a greater depth. Propagation of 
this type is shown in Figure 13. 

3) The angle between the zero degree ray and the 


first ray to touch bottom (trapping angle) was greater in 
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the cold-to-warm direction in 18 out of the 21 pairs of ray 
traces which were computed. The exceptions were the 3 
pairs of runs from Section E which had the smallest surface 
thermal gradient of any of the sections. 

4) A large percentage of rays propagated in the cold- 
to-warm direction attain their maximum depth at or near the 
maximum range of the run. This indicates that the cold-to 
Warm rays are diving deeper and deeper with increasing range. 

5) As the surface thermal gradient is increased, ray 
traces in opposing directions become more similar in both 
surface bounces and ray depth. Tnis is primarily a result 
of the interpolation of sound velocity from one profile to 
another, 

Of particular interest are the curves of Figure 16. 
The surface thermal gradient is plotted versus the differ- 
ence in angle trapped for each pair of runs in opposite 
directions, The convention used is that the angle obtained 
by subtracting the trapping angle of the warm-to-cold run 
from the trapping angle of the cold-to-warm run is con- 
Sidered to be positive. These curves, drawn for each of 
the assumed profile spacings, indicate that the difference 
in angle trapped tends to increase as the surface gradient 
increases for each of the three distances examined. At 
first glance, the decrease in the difference as the spacing 
of the sections decreases from 10 to 6 to 4 xp appears to 


be in contradiction to the basic premise that the difference 
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in trapping angle increases as the surface gradient in- 
creases. However, close examination of the angles trapped 
shows that as the profiles are moved more closely together, 
the trapping angle in the warm-to-cold direction of propa- 
gation increases. [It then more closely approaches the cold- 
to-warm trapping angle which remains relatively stable. 
This occurs in six of the seven sections examined. Only. 
Section C does not exhibit this tendency at a profile 
spacing of 4 KM. It has not been determined why this pair 
of runs is inconsistent with the remainder of the data. 
Thus, Figure 16 points to the possibility of preparing 
operationally useful nomograms for estimating sonar 
efficiency in frontal regions. These nomograms would 
require the collection of more closely spaced synoptic 
observations than were available to this author but could 


be of immense operational use in oceanic frontal zones. 
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V. CONCLUSIONS 

Computation of sonar ray traces across the oceanic 
polar front in the Denmark Strait region and along the 
southeastern coast of Greenland and an examination of the 
hydrography of the area lead to the following conclusions; 

1) The position of the oceanic polar front may vary 
Widely from day to day, but generally is found along the 
edge of the continental shelf. 

2) While local influences may temporarily alter the 
gradient structure, conditions on an annual basis result 
from the effect of atmospheric wind systems which drive 
the Irminger and East Greenland Currents, 

3) The angle of propagation trapped by sonar rays 
propagating across the front is larger when the sound 
source is in cold water, compared to conditions with the 
source in warm Water. 

4) As the gradient increases, the angle trapped by 
propagation from the Warm side increases. Gradient has 
little effect on the cold-to-warm trapping which varies 
only slightly. 

5) A relationship can be demonstrated between the 
surface thermal gradient and the change in angle trapped 
with change in direction of propagation for each assumed 


Preoerile spacing. 
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VI. RECOMMENDATIONS FOR FURTHER RESEARCH 

A requirement exists for a thorough examination of 
this frontal zone by towed thermistor chain or closely 
Spaced XBiepror miles. 

Repetitive sections should be taken in the same lo- 
cations as often as possible, but at least once during each 
season. This would provide an excellent data base for the 
conetruction of the sonar efficiency nomograms mentioned 
in Seerion iy; 

An additional objective should be clarification of the 
relationship of the strength and position of the oceanic 


polar front to the atmospheric polar front. 
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